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A FOURTH-ORDER MODEL FOR MEMS WITH 
CLAMPED BOUNDARY CONDITIONS 

PHILIPPE LAURENgOT AND CHRISTOPH WALKER 



(^ I Abstract. The dynamical and stationary behaviors of a fourth-order equation in the unit ball 

CN ■ with clamped boundary conditions and a singular reaction term are investigated. The equation 

$H i arises in the modeling of microelectromechanical systems (MEMS) and includes a positive volt- 

^H. age parameter A. It is shown that there is a threshold value A* > of the voltage parameter such 

that no radially symmetric stationary solution exists for A > A*, while at least two such solutions 

QQ ■ exist for A G (0, A*). Local and global well-posedness results are obtained for the corresponding 

hyperbolic and parabolic evolution problems as well as the occurrence of finite time singularities 

when A > A*. 



< 



1. Introduction 

> ■ 
\Q , Electrostatically actuated microelectromechanical systems (MEMS) are microscopic devices 

Gn I which combine mechanical and electrostatic effects. A typical MEMS device is made of a rigid 
i!^ ' conducting ground plate above which a clamped deformable membrane coated with a thin con- 
ducting film is suspended. Application of a voltage difference induces a Coulomb force which, 
(^ '. in tum, generates a displacement of the membrane. An ubiquitous feature of such devices is, 
(^ I that when the applied voltage exceeds a certain threshold value, the membrane might collapse 
(or touch down) on the ground plate. Controlling the occurrence of this phenomenon - usually 
referred to as the "pull-in" instability - is of utmost practical importance in the design of such 
^ . devices either to set up optimal operating conditions or to avoid device damaging. Mathemati- 
j^ I cal models have been derived to describe MEMS devices which lead to free boundary problems 
due to the deformable membrane [|2n . Since these models are difficult to analyze mathemati- 
cally (though recent contributions can be found in [|7l HI [151), one often takes advantage of the 
small aspect ratio of the devices to reduce the free boundary problem to a single equation for the 
displacement, see lETl . More precisely, the small aspect ratio model describing the dynamics 
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of the displacement u = u{t, x) of the membrane f2 C M'^ reads 

-f^d^u + dtu + BA^u - TAu 
u = duU 

Here, "f'^d'^u and dtU account, respectively, for inertia and damping effects, BA'^u and —TAu 
are due to bending and stretching of the membrane, while — A(l + m)~^ reflects the action of 
the electrostatic forces in the small aspect ratio limit. The parameter A is proportional to the 
square of the applied voltage. Observe that the right-hand side of (11.11) features a singularity 
when n = — 1, which corresponds to the touchdown phenomenon already mentioned above. 
Since the strength of the singular reaction term is tuned by the parameter A, it is not surprising 
that the latter governs the existence of stationary solutions, that is, solutions to 

BA\-TAu = -- -, xen, (1.4) 

(1 + u)^ 

u = d^u = 0, xedQ. (1.5) 

When bending is neglected, that is, when B = 0, this problem reduces to a second-order elliptic 
equation that has been studied extensively in the recent past, see e.g. the monograph [ITOl and 
the references therein. As expected from the physics there is a critical value A^, > such that 
no stationary solution exists if A > A* and at least one stationary solution exists for A G (0, A*). 
Let us emphasize that the comparison principle is available in this case and turns out to be a key 
tool for the analysis. Less attention has been dedicated to (ll.4l) - (ll.5l) with B > 0, one reason 
might be the lack of a maximum principle in general for the clamped boundary conditions (11.51) 
(also called Dirichlet boundary conditions), see the monograph [[TT| for a detailed discussion 
of positivity properties of higher-order operators. Recall that the situation is completely dif- 
ferent if the clamped boundary conditions (11.21 ) are replaced with pinned (or Navier) boundary 
conditions 

u = Au = on on, (1.6) 

since in this case the maximum principle holds in arbitrary domains [fTOlfTTIl . This allows one 
in particular to show similar results for the fourth-order problem (11.41) . (11.61) as outlined above 
for the second-order case corresponding to B = 0. We refer to [[TOl [TTl for details. 

Returning to the case of clamped boundary conditions, when B > existence of solutions 
to (ll.4l) - (ll.5l) for small values of A has been established in ifTTl for an arbitrary domain il. This 
is the only result we are aware of for a general domain. In the particular case when fi equals 
the unit ball Bi, Boggio ^ has uncovered the availability of the maximum principle for the 
operator BA"^ with boundary conditions (11.51 ) by showing that the corresponding Green function 
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is positive. This fact has been used in ifTOl to describe more precisely the set of solutions to (11.41 )- 
(11.51) with T = 0. It has actually been shown in ifTOl Chapter 11] that there is a critical threshold 
value A* > such that no solution exists for A > A* and a solution exists for A G (0, A*]. 

Very recently we were able to extend Boggio's maximum principle to the operator BA^ — TA 
with T > and boundary conditions (11.51 ) in the class of radially symmetric functions in Bi 
|[T6l . Taking advantage of this property not only allows us to extend the results of lITOl to include 
T > (for d = 1,2), but also to show that for each voltage value A G (0, A*) there are at least 
two (radially symmetric) solutions to (ll.4l )- (ll.5l ). thereby answering a question raised in [ITOl 
p. 268]. A summary of our results on radially symmetric solutions to (ll.4l) - (ll.5l) is stated in the 
following theorem. 

Theorem 1.1. Let d = 1,2, Q = Mi, B > 0, and T > 0. There exists A^, > such that there 
is no radially symmetric solution to (|1.4| )- (|1.5I) for A > A^.. Moreover, there is a continuous 
curve (A(s), [/(s)), s G [0, oo) in M x C^(]Bi) such that U{s) is for each s G [0, oo) a ra- 
dially symmetric solution to (ll.4l) - (ll.5l) with A = A(s). Moreover, (A(0),f/(0)) = (0,0) and 
(A(s), U{s)) — )■ (0,0;) as s -^ oo, where u is an explicitly given radially symmetric function 
with Ci;(0) = —1. Finally, there is s* > such that A is an increasing function from [0, s*] onto 
[0, A^,] and A is decreasing in a right-neighborhood of s^.. 

Remark 1.2. Theorem 17.71 guarantees that for each A G (0, A^,), there are at least two (radially 
symmetric) solutions to (|1.4| )- (|1.5I ). More precisely, for each A G (0, A^,) there are at least 
two values < Si < s^, < S2 with A(sj) = Xforj = 1,2 and U{s2) < U{si) in Bi with 
U{S2) ^ Uisi). 

Let us mention here that the construction of the curve (A(s), U{s)), s G [0, s*] follows the 
lines of lUOl Chapter 11], where a similar result is proved when T = 0. There are, however, 
some technical difficulties to be overcome. Nevertheless, we emphasize that the main contribu- 
tions of Theorem 1 1 . 1 1 are the extension of the curve past (A(s*), f/(s*)) and the identification of 
its end point w as s — )■ oo. Interestingly, the end point to is given as a solution of a boundary 
value problem in Bi \ {0} which can be computed explicitly (see Theorem 12. 20| below); a plot 
of u is shown in Figure [T] The qualitative behavior of cu is the same for d = 1 and d = 2. 

For the case of pinned boundary conditions (11.61 ) it has been shown in |fT3l Theorem 1.2] with 
the help of the Mountain Pass Principle that there are at least two solutions for A G (0, A*). The 
limit as A —7- of the minimum of the solutions constructed with the Mountain Pass Principle 
is proved to be —1. However, the precise profile as A —> is not identified therein. 

The proof of Theorem II. H is performed in Section |2] Therein we give a more detailed char- 
acterization of the set of radially symmetric stationary solutions. Actually, the implicit function 
theorem provides a branch ^o of radially symmetric solutions (A, u) to (|1.41 )- (|1.51 ) emanating 
from (0, 0). We then use the bifurcation theory of [|5l for real analytic functions to extend ^o to 
a global curve A (see Theorem |2.5| below). Next we show that ^o coincides with the branch of 
stable radially symmetric stationary solutions (see Corollary 12.161) . To achieve this result, the 
maximum principle obtained in liT6l is essential as was Boggio's maximum principle in [JOj for 



Ph. LAURENgOT AND CH. WALKER 




Figure 1 . Plot of the function w for d = 2 and (S, T) = (1, 50). 

the case T = 0. The outcome of this analysis is that there is a threshold value A* > such that 
there is no radially symmetric stationary solution for A > A*, while for any A G (0, A*) there is a 
unique stationary solution ux such that (A, ux) E Aq. We then show that for A = A^, there is also 
a radially symmetric stationary solution ux, , which guarantees on the one hand that Aq ^ A 
and on the other hand that we may apply the result of DHl to extend the branch ^o to the "right" 
of (A*, maJ (see Theorem [2.18l) . The final step is to show that A connects (A, u) = (0, 0) to the 
end point (0, u) and to identify the latter (see Theorem 12. 201 ) . As a consequence, the continuous 
curve A passes through (A, u) = (A*, maJ, which implies Remark [L2l 



We shall also investigate local and global well-posedness of the dynamic problem (ll.lj l- 
(11.31) . It is worth pointing out that the maximum principle, which is at the heart of the proof of 
Theorem II. II and the main reason to restrict the analysis to Bi, is no longer valid for the time- 
dependent problem (ll.ll) - (ll.3l) . In order to construct solutions to the evolution problem, we 
therefore have to employ an alternative method which does not rely on maximum principles. 
Our approach is based on semigroup theory and is not specific to Bi. We thus consider an 
arbitrary domain Q C W^, d = 1,2, in the following and begin with the hyperbolic problem 
which has not received much attention so far. 



Theorem 1.3. Let Q G M.'^ be an arbitrary bounded smooth domain for d = 1,2 and B > 0, 
T>0. Let\>OandKe (0,1). Let{u^,u^) e H^{Vt) xH^{Vt) be such that u^ >-l + Kinn 
and such that n° and u^ both satisfy the boundary conditions (|1.2I) . Then the following hold: 

(i) There are Tm > and a unique maximal solution u to (ll.ll) - (ll.3l) with regularity 

n G C([0, rj, i72(fi)) n ^^([0, rj, L2(fi)) 
for /c = 0, 1, 2 and r G (0, r^). 



d1ueL^{Q,T;m-'\^)) 
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(ii) If Tm < oo, then 

liminf (minM(t)) = —1 . 

(Hi) There are Ai(k) > and r^n) > such that Tm = oo provided that A < Ai(k) and 
\\{u^,u^)\\jj2y^x.2 — ''^{i^)- In Ms case, u G Loo(0, oo; H^{Q)) and 

inf M > — 1 . 

(O,oo)xn 

(iv) IfQ = Bi and both u^ and u^ are radially symmetric, then so is u{t)for each t G [0, Tm). 

Similar results have been established in [|T4l for rf = 1 and B = (without the damping term 
dfU) and in lfT2ll for the pinned boundary conditions (11.61 ) and d E {1, 2, 3}. Let us point out 
that the semigroup approach allows us to obtain strong solutions instead of weak solutions as in 
HH m, see Section O for the proof of Theorem [O] 

In the damping dominated limit 7^1 when viscous forces dominate over inertial forces, 
(ll.ll) - (ll.3l) reduces to the parabolic problem 

dtu + BA'^u - TAu = 



>0, 


X e^, 


(1.7) 


>0, 


X edn, 


(1.8) 




X eQ . 


(1.9) 



[i + uy ' 

u = dyU = , 

u{0,-) = u\ 

To the best of our knowledge, this problem has not been investigated so far. For this reason, we 
include a result on its well-posedness though local existence of solutions is a rather classical 
argument. To obtain global solutions for small values of A, we consider only regular initial 
values in the next theorem for the sake of simplicity. 

Theorem 1.4. Let Q G M.'^ be an arbitrary bounded smooth domain for d = 1,2 and B > 0, 

T > 0. Let X > and k e (0, 1). Let n° G H'^{^) be such that n° > -1 + k in Q and vP 
satisfies the boundary conditions (11.81 ). Then the following hold: 

(i) There are r^ > and a unique maximal solution u to (ll.7l) - (|1.9l ) with regularity 

u G c([o, Tm), H\9)) n c((o, Tm), H^Q)) n c\{o, Tm), L2{n)) . 

(ii) If Tm < oo, then 

liminf (minu(t)) = —1 . 

(iii) There is \\{^k) > such that Tm = oo provided that A < Ai(k) and \\u^\\fj2 < k~^. In 
this case, u G Loo{0, oo; H'^(Q)) and 

inf M > — 1 . 

(0,oo)xr2 

(iv) If Q = Ml and u^ is radially symmetric, then so is u{t) for each t G [0, Tm). 
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The proof of Theorem II .41 follows the same lines as that of Theorem 11.31 and is to be found 
in Section [3]2l 

On physical grounds it is expected that solutions to the dynamic problems (ll.ll) - (ll.3l) or (11.71 )- 
(11.91 ) touch down (i.e. u = —1) in finite time and thus do not exist globally if the voltage value 
A exceeds the critical pull-in voltage above which no stationary solution exists. This is true for 
the second-order parabolic case, see [jTOl for instance, but seems to be an open problem both for 
the hyperbolic equation (ll.ll) - (ll.3l) as well as for the parabolic equation (ll.7l) - (|1.9l ). Actually, 
even the weaker result of the occurrence of touchdown in finite time for large values of A has 
not yet been proven apparently, though observed numerically in [[TSl for (ll.7l) - (|1.9l ) with T = 
and shown in [|T4]| in the absence of bending {B = 0). The next result is a step in that direction 
when Q is the unit ball Bi of R'^. 

Proposition 1.5. Let 1] = Bi with d e {1,2}, X > 0, B > 0, T > 0, and let (n°, n^) G 
if^(Bi) X if ^(Bi) be such that vP > —1 in Bi and both u^ and u^ satisfy the boundary con- 
ditions (11.21 ). Let u be the maximal solution to either (|1.1| )- (|1.3I ) with initial condition {vP, u^) 
enjoying the properties listed in Theorem \L3\ or (|1.7| )- (|1.9I) with initial condition n° enjoying 
the properties listed in Theorem \L4\ Let t„i be its maximal existence time. If A is sufficiently 
large (see (14.21) below for a quantitative lower bound), then t„i < oo. 

It is worth pointing out that the outcome of Proposition 11.51 complies with the numerical 
simulations of (ll.7l) - (|1.9l) performed in lITSl in Bi and showing the occurrence of finite time 
touchdown. The proof of Proposition II .5 1 is given in Section l4~n and relies on the eigenfunction 
method. 

Owing to the study carried out in Section |2l we are able to refine this result in the radi- 
ally symmetric setting and show that the touchdown behavior indeed starts exactly above the 
threshold value A* defined in Theorem ll.il 

Proposition 1.6. Assume fi = Bi with d G {1,2} and let (m°, u^) be radially symmetric initial 
conditions satisfying the requirements of Theorem \L3\ ifi > or Theorem \L4\ if i = 0. Then, 
if X > A*, the corresponding maximal solution to (|1.1I) - (|1.3I) or (|1.7I) - (|1.9I) on [0, Tm) does not 
exist globally, that is, Tm < oo. 

The proof of Proposition ll .6l is performed in Section 143] and also relies on the eigenfunction 
method, but with a more accurate choice than in the proof of Proposition 1 1.5 1 as already noticed 

in ma. 

Let us conclude the introduction with some remarks on the qualitative behavior of solutions to 
the evolution problem in the ball Bi. Proposition 1 1.5 1 and Proposition II .61 show the occurrence 
of a finite time singularity, but do not provide information about the precise behavior near 
touchdown time. According to the numerical simulations performed in lITSl . the fourth-order 
term has a strong influence on the way solutions touch down in finite time as this might take 
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place on a circle (for d = 2). This markedly contrasts the second-order case, where touchdown 
occurs only at the single point a; = 0, see ifTOl Theorem 8.3.4]. 

When a solution to (ll.ll) - (ll.3l) does not touch down in finite time, then it exists globally in 
time and might even be bounded away from —1 as well as be bounded in H^ according to 
Theorem 11.31 (if 7 > 0) and Theorem 11.41 (if 7 = 0). A natural next step to understand its 
dynamics is to investigate its large time behavior. While this seems to be an open problem for 
a general domain ^2, the analysis performed in Section |2]for fi = Bi in the radially symmetric 
setting paves the way for a better understanding of this issue. On the one hand, Proposition l2.15l 
below entails that one may apply the principle of linearized stability to show that U{s) is locally 
asymptotically stable when s G [0, s*). On the other hand, it might be possible to use the 
Lojasiewicz-Simon inequality as in ||T21 to establish convergence to a single steady-state. 

2. Radially symmetric stationary solutions 

In the following, if S{i^) is a space of functions defined on f2, we write Soi^) for the sub- 
space of S{Q) consisting of functions u satisfying the Dirichlet boundary conditions (11.51 ). if 
meaningful. If f2 is the unit ball Bi of M"', then S'r(Bi) stands for the subspace of S'(Bi) con- 
sisting of radially symmetric functions. Clearly, SD,ri9i) '■= Sd{^i) H S'r.(Bi). If S{Q,) is a 
normed vector space, then || ■ \\s stands for its norm. For p G [1, 00] we denote the norm of 
Lp{il) simply by || ■ ||p. 

Recall that the stationary solutions of (|1.1| )- (|1.2| ) satisfy 

BA'^u-TAu = -Xg{u) in Bi , (2.1) 

u = d^u = on aBi , (2.2) 

where g{^) := (1 + 0"^ for ^ > -1. 

Definition 2.1. A radially symmetric classical solution u (with parameter X) of the boundary 
value problem (|2.1| )- (|2.2| ) is a radially symmetric function u G C^(Bi) fl C^(Bi) satisfying 
u > —1 in Bi and solving (|2.1| )- (|2.2| ) in the classical sense. 

We denote the set of all radially symmetric classical solutions with parameter A to the bound- 
ary value problem (|2.1I) - (|2.2I) by S^. 

Similarly, a radially symmetric function u G C^(Bi) n C^(Bi) satisfying ?i > — 1 in Bi is a 
classical subsolution of (I2.1I) - (I2.2I) (with parameter A), if it satisfies B/S?u — TAu < —\g{u) 
inBianddlJlonSBi. 

We introduce the operator 

Au := BA^u - TAu , u G i^D,r(®i) , 
and recall the following well-known properties: 
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Lemma 2.2. A e £(if|, ,.(Bi), L2,r(Bi)) is invertible with 

for each a G (0, 1). Moreover, there are mi > and 0i G Cj^^i^i) with 0i > m Bi, 

||0i||i = 1, and A(j)i = mi0i. 

Proof. The invertibility of A and the regularity properties of A^^ are consequences of [fTTl 
Theorem 2.15, Theorem 2.19, Theorem 2.20]. That there is a positive eigenvalue with a positive 
eigenfunction follows from [[T6l Theorem 4.7]. D 

We further define 

A* := sup {A > : 5^ is non-empty} G [0, oo] , (2.3) 

and first derive some elementary properties of S^. 

Lemma 2.3. The following hold: 

(i) 5° = {0} andS^ C C|.,(ii) /or A > 0; 
(ii) ifX > and u G S^, then —l<u<Oin Bi; 
(iii) the threshold value A* defined in (12. 3|) is finite. 

Proof. The first statement of (i) readily follows from Lemma |Z21 If m G S^, then (?(«) belongs 
to C^(Bi) since (? is smooth in (—1, oo) and m > — 1 in Bi. Thus u G C|,^(Bi) by Lemma [Z2l 
Moreover, —1 < u < in Bi by 11161 Theorem 1.4] since \g{u) > (see also Lemma 1277] 
below). Consequently, testing (12. II) - (12. 21) by 0i > introduced in Lemma IZ2] yields 

—mi / (f)idx < nil (f)iudx = —A / (pi g{u) dx < —A / 01 dx , 

JBi J^i Jbi Jbi 

whence A < mi. Therefore, A* < mi < oo. D 

In fact, one can show that A* < mi. Indeed, assume A* = mi for contradiction so that there 
are sequences A„ — )■ mi and m„ G S^". Then, the above computation actually yields 



Since also 



mi- Xn> A„ / (pi {g{un) - 1) dx > 



(t)i{g{un)-l)dx= / 01 " " ,„ dx > / 0i|u„|dx 

(1 + Mn)^ 



we conclude that 



lim A„ / (pi g{un) dx = nil and lim / 0i |m„| dx = . 

'^^'^ JBi "-^"^ ill 



A FOURTH-ORDER MEMS MODEL WITH CLAMPED BOUNDARY CONDITIONS 9 

From Un E S^" we obtain 

mi (t)iUndx = -Xn (j)ig{un)dx, 

and letting n — )• oo and using the previous limits imply that mi = 0, which contradicts 
Lemma [X2I 

Remark 2.4. Observe that the computation in the proof of Lemma 12.31 excludes also the exis- 
tence of non-radially symmetric solutions to (II ■4| )- (|1 .51 ) for A > rrti. An interesting question is 
whether there are non-radially symmetric solutions for A G (A*, mi). This is not the case when 
T = as it is shown in [|2l that all solutions to (ll.4l) - (ll.5l) are radially symmetric. 

2.1. A continuous curve of stationary solutions. In this subsection we invoke the global bi- 
furcation theory of [5, Section 2.1] for real analytic functions to establish the existence of a 
global curve of radially symmetric stationary solutions. This tool has also been used in [|TOl 
Section 6.2] for the second-order case (that is, B = 0). 
Since 

) : |m| < linii} (2.4) 
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F: 


RxO^CJj, 


.(Bl) 



{\,u)\^u + \A~^g{u) , (2.5) 

is well-defined according to Lemma lZ2] and real analytic. Observe that n G S^ if (A, n) G Mx (9 
with F(A, u) = 0, the bound |n| < 1 following from Lemma |Z3] Clearly, F(0, 0) = and the 
partial Frechet derivative Fu(0, 0) equals theidentity inC]5^(]Bi). Thus, by the implicit function 
theorem, the zeros of F near (0, 0) are given by a real analytic curve (A, V{X)) with V^(0) = 0. 
Moreover, there exists Aq G (0, cxd], which is maximal with respect to the existence of a real 
analytic function V : [0, Aq) ^ C}j^^{Mi) for which F(A,1/(A)) = and i^,,(A, 1/(A)) G 
C{C}j^^{Mi)) is boundedly invertible for each A G [0, Aq). Consequently, the set 

S := {(A, u) G (0, 00) X C : F(A, n) = and F„(A, u) G £(C^ ,,(Bi)) 

is boundedly invertible} 

is non-empty as it contains the maximal arc-connected subset 

Ao:={{\,V{X)) : Ag(0,Ao)}. (2.6) 

Note that Aq and V are unique and necessarily Aq is finite since it belongs to (0, A^,]. We have 
thus verified assumption (CI) from dH Section 2.1]. For (C2) therein we may argue as in [[TOl 
p. 128] that this assumption merely serves to show in the proof of [|51 Theorem 2.3 (iii)] that 
5* is open in its closure S and can thus be replaced by the stronger one that (0, 00) x C is 
open in M X C}j ^.(IBi). Then, since Hj, ,.(Bi) embeds compactly in C}) ^(IBi), we may regard 
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the operator XA^^g'(u) E C{C\) ^(Bi), Hj^ ,,.(Bi)) as a compact operator in C\) r(IBi) for each 
(A, m) G M X O. Hence, by ll22l Theorem 4.25], the partial Frechet derivative 

Fu{\ u) = l + \A-^g'{u) , (A, m) e M X O , 

is a Fredholm operator of index 0. The remark in []5] p. 246] now entails that (C3)-(C5) therein 
hold. Next, we introduce the function 

i^ : (0, oo) X C -)■ [0, cx)) , (A, u) t-). 



min^Jl + n} ' 

To verify (C6) from [JS] consider a sequence (A„, Un)neN in (0, oo) x O with F(A„, m„) = and 
i/(A„, Un) < c < oo for each n G N. Then, by Lemma [Z2l each m„ belongs to C|, ,.(lBi) with 
Un> —I + c~^ in Bi and satisfies 

BA\n - TAun = -Xn9{un) in Bi . 

The above uniform lower bound on n„ and the finiteness of A* established in Lemma [Z3l now 
imply that the sequence {Xng{un))neN is bounded in Loo(Bi). Thus, {Xn,Un)neN is bounded 
in (0, A^,] X Hf) ^(Bi) and so (A„, Un)neN has a converging subsequence in [0, A^.] x C}j ^(Bi). 
Hence (C6) in Q holds true. Setting A := 0, we clearly have (A, 0) ^ (0, oo) x C and (A, 0) 
is in the closure of ^o defined in (12.61) . whence (C7) in ^. Finally, suppose that (A„, Mn)nGN 
is a sequence in (0, oo) x O with F{\n,Un) = and z/(A„,m„) < c < oo for each n E N, 
which converges in M x Cjj ^(^i) towards (A, u) ^ (0, oo) x O. Then — 1 + c~^ < Un < in 
Bi and < A„ < A* for each n G N by Lemma |23l which entails that — 1 + c~^ < u < 0, 
whence (A, u) G [0, A*] x O and F(A, u) = 0. Since (A, u) ^ (0, oo) x O, this is only possible 
if (A,m) = (0,0). The implicit function theorem guarantees that (A„,m„) G ^o for '^ large 
enough. This yields (C8) in [|51, and therefore, we are in a position to apply [|5l Theorem 2.4] 
and obtain: 

Theorem 2.5. There is a continuous function (A, U) : (0, oo) — )■ (0, oo) x C})j.(Ei) with the 
following properties: 

(i) U{s) G Sr^'^ for each s G (0, oo); 
(ii) (A,t/)((0,1)) C Aa«^lim,_o(A(s),t/(s)) = (0,0); 
(iii) (A, U) is injective on (A, U)~^{S) and 

lim (niinf/(s)) = — 1 ; 

(iv) at all points s G (A, U)~^{S), (A, U) is real analytic with A'(s) ^ 0. 
Actually, more precise information is given in [|5] Theorem 2.4] about the curve 

A:={{K{s),U{s)) : sG(0,oo)}, (2.7) 

traced out by the function (A, f/), in particular, that it is piecewise analytic: 
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Remark 2.6. The set (A, U)~^{S \ S) C. (0, oo) consists of isolated values and locally near 
each point sq G (A, U)~^{S \ S), there is a re-parametrization ( of the parameter s such that 
(A, U) o ( is real analytic with derivative vanishing possibly only at 0. 

Before analyzing further the curve A and in particular showing that it "globally" extends ^o 
defined in (12.61) (note that at this point, the curves ^o and A could still coincide), we first derive 
general properties of solutions to (|1.41 )- (|1.51 l in the next subsection . 

2.2. Properties of stationary solutions. We first recall the following sign-preserving property 
of the operator BA^ — TA with homogeneous clamped boundary conditions in Bi with radial 
symmetry established in [[T6l . 

Lemma 2.7. Consider/ G Cr{Mi) andw E C^ (Mi) nCf. (Mi) such that w is a classical solution 
to 





BA^w-TAw = f in Bi , 




w = dyW = on (9Bi 


Then, if f < in Mi, 




Furthermore, 


either w = or w < in B 



(2.8) 
minw = w(0) , (2.9) 

Bi 

and there is vq G (0,1) such that 

Aw<0 in Bi\B^„ and Aw > in B,^ \ {0} . (2.10) 

Finally, the profile ro ofw defined by tt>(|a;|) = w(x) for x G Bi is a non-decreasing function 

on [0, 1]. 

Proof. The first statement (12.81 ) of Lemma 12771 readily follows from [[161 Theorem 1.4]. Fur- 
thermore, the proof of that result reveals that ([2. lot is true. We next deduce from (12.101) 
that dr (r'^''^drK){r)) < for r G (ro, 1] and dr (r^~^drKi{r)^ > for r G (0,ro). Since 
drtv{0) = 9rlt)(l) = due to the radial symmetry of w, its regularity, and its boundary condi- 
tions, we conclude that drW{r) > for r G [0, 1]. Then ro is a non-decreasing function in [0, 1] 
and attains its minimum at r = 0. D 

Lemma 2.8. Define the scalar product (■, ■) on H^) ^(Bi) by 

{v, w) := / [BAv{x)Aw{x) + TVv{x) ■ Vw{x)] dx , v,w e Hl^{Mi) . 

Let /C := {v E H^) ^(Bi) : v > O} be the positive cone of Hjj ^(Bi) and define its polar cone 
by 

/C° := {w E HlX^i) : {v,w) < for all v E /C} . 

Then, given v E H^y ^(^i)' there is a unique couple (fi, V2) G /C x /C° such that {vi, V2) = 
and V = vi + V2- In addition, V2 < a.e. in Bi. 
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Proof. The fact that any v E iJ|, ^(Bi) can be written in a unique way as a sum f = wi + t;2 
with {vi, V2) = and {vi, t>2) G /C x /C° is a well-known result due to Moreau [|T9l . The non- 
negativity property of V2 actually follows from the sign-preserving property stated in Lemma [277] 
and can be proved as [|T6l Proposition 4.5], where only the one-dimensional case was handled. 

D 



As in [flOl the linear stability of stationary solutions is an important tool in the detailed anal- 
ysis to follow. For M G iS^, it is measured by 

/ii(M) := inf I f {B\Av\^ + T\\/v\^ + Xg^u^) dx : v E Hl^^(Ei) , \\v\\2 = 1 j , 

(2.11) 
which turns out to be a simple eigenvalue of the linearization of (|2.1I) when non-negative as 
shown in the following lemma. 

Lemma 2.9. Consider A G [0, A*] and u E S^ such that Hiiu) > 0. Then the following hold: 
(i) fii{u) is a simple eigenvalue of the operator A + \g'{u) G C{Hf) ri^i), -^2,r(Ki)) with 

a positive eigenfunction in C|, ^(^i)'' 
(ii) fii{u) > if and only if Fu{X,u) = 1 + XA^^g'{u) G C{C}j ,^{Mi)) is boundedly 
invertible. 

Proof, (i) A classical compactness argument along with the weak lower semicontinuity of the 
scalar product (■, ■) in H^ r(^i) defined in Lemma lZSl guarantee the existence of a minimizer 
to (12.111) in Hl)^.{Mi) satisfying ||0||2 = 1. Then G H^^{Mi) is a solution to the correspond- 
ing Euler-Lagrange equation 

fiA20-TA0+(A^'(n) -/ii(n))0 = O in Bi , = 9^0 = on ^Bi . (2.12) 

Now, let G Hfjj.(Mi) be any solution to the boundary value problem (12.121 ). According to 
Lemma [Z8l there is a unique couple (0i, ^2) G /C x /C° such that = 0i + 02. (0i7 02) = 0, and 
02 < a.e. in Bi. We deduce from the definition (|2.11l) of yUi(n), the orthogonality properties 
of (01, 02), and (EJIl) that 



/il(M)||01 -02II2 < (01 -02,01 -02) + A / 5''(m)(01 -02)^dx 

< (01 + 02, 01 + 02) + A / g'{u){4>i - 4>2f dx 

< A / g'{u) [(01 - 02)=^ - (01 + 02)'] dx + /ii( 



u 



^2||2 



whence 



< — 4A / g'{u)(f)i(f)2 dx + 4/ii(M) / 0i02 dx . 
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Both terms of the right-hand side of the above inequality being non-positive, we infer from the 
negativity of g' that 

0102 = a.e. in Bi . (2.13) 

Now, for i = 1, 2, it follows from the embedding of H^{Mi) in C"(Bi) for a e (0, 1) (recall 
that d e {1, 2}) that 0j G C"(Bi) and, according to [|TT1 Theorem 2.19], the boundary value 
problem 

BA'^ipi - TAtpi = [iJi{u) - Xg'{u)] (pi in Bi , V'i = duipi = on ^Bi , (2.14) 

has a unique classical radially symmetric solution ipi G Cj^"*'"(Bi). Since /ii(n) — \g'{u) > 0, 

01 > and 02 < in Bi, it follows from LemmaOthat V'l > > ^2 in Bi with ^1 > in Bi 
if 01 ^ and ^2 < in Bi if 02 ^ 0. In addition, due to dTHl) and (I04l) . 

5A2(0-^i-V;2)-TA(0-^i-^2) = K(n)-A(?'(n)](0-0i-02) = O in Bi 

with — ^1 — ^"2 £ Hf)^(Mi), whence = ^1 + ip2. Furthermore, ipi clearly belongs to /C 
while, for any w G /C, we infer from (12.141 ) that 

(^^2, V) = ifJ'liu) - Xg'{u)) 0217 dx < , 

SO that 02 £ ^°- The uniqueness of Moreau's decomposition then warrants that ^j = 0j for 
i = 1,2. Therefore, if 0i ^ and 02 ^ 0, we deduce from the above analysis that ipiip2 < 
a.e. in Bi and ^1^2 = 0i02 = a.e. in Bi, and a contradiction. Therefore, either 0i = or 

02 = 0, and we have shown that does not change sign in Bi. 

Consequently, any element of the kernel of the operator A + \g'{u) — yUi(u) in H'^^iMi) 
does not change sign, which implies that the kernel's dimension is one by a classical argument. 
Indeed, assume for contradiction that there are two linearly independent positive functions 
and ip in the kernel. Then — aip with suitable a > is a sign-changing function in the kernel, 
which is impossible. Therefore, the kernel of A + \g'{u) is spanned by a positive function 
G Cf) r(®i)' ths additional regularity stemming from Lemma |2]2] Finally, to show that /ii(7i) 
is a simple eigenvalue oi A + Xg'{u), consider $ G i^l, ^(Bi) such that A^ G Hfj^.{Mi) and 
{A + Xg'{u) - i2i{u))'^^ = 0. Then, {A + Xg'{u) - fii{u))^ = a0 for some a G M. Multiplying 
this identity by and integrating over Bi gives a||0||| = 0, thus a = 0. This yields assertion (i). 

(ii) Assume that F„(A,n) = 1 + \A'^g'{u) G £(C|)^(Bi)) is not boundedly invertible. 
Then —1 is an eigenvalue of the compact operator \A'^g'{u) G C{C}j .^(Mi)). Hence there is 
G Ci,^(Bi) with + AA-i^'(m)0 = 0. Altematively, A0 = -A^'(w)0 so that G Hj)^^{Mi) 
by Lemma [2!2l and 

/ii(M)||0||^<(0,0) + A / (7'(m)0Mx = O, 



which implies /ii(n) < 0. Conversely, if /ii(n) = 0, then, arguing as in the proof of Lemma 12. 91 
we obtain a solution G Hfj^(Bi) to 

BA'^cf) - TA0 + A^'(n)0 = in Bi , 
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and thus (j) + XA~'^g'{u)(j) = 0. D 

As in lUOl Chapter 1 1], a key tool in the analysis is the following comparison lemma. 

Lemma 2.10. Consider A G (0, A*] andu G S^ such that niiu) > 0. 

(i) Ifve C|) ,,(Bi) n C^(ii) is a classical subsolution to (l2?T1)-(IT2l) with v > -I in li, 

then V <u in Bi. 
(ii) Furthermore, v = u if Hiiu) = 0. 

Proof, (i) We proceed along the lines of the proof of [[TOl Lemma 11. 3. 4]. By Lemma IT8] there 
is a unique couple {wi,W2) G /C x /C° such that v — u = Wi + W2, (wi, W2) = 0, and W2 < 
a.e. in Bi. Since 

BA'^iv -u)- TA{v -u) + X{9{v) - g{u)) < in Bi 

and wi G /C, we may multiply the above inequality by wi and integrate over Bi to obtain 

{wi, V — u) + \ / {g{v) — g{u))wi dx < . 

ill 

We next deduce from /ii (n) > that 

{wi, V — u) = {wi, Wi) > —A / g'{u)w\ dx = —A / g'{u)wi{v — u — W2) dx 

= ~^ / g'{u)wi{v — n) dx + A / g'{u)wiW2 dx . 
Combining the previous two inequalities gives 

A / {g{v) — g{u) — g'{u){v — u))wi dx + X / g'{u)wiW2dx < . 



Owing to the convexity and the monotonicity of g together with the sign properties of wi and 
W2, the two terms on the left-hand side of the above inequality are non-negative. Therefore, 

{g{v) — g{u) — g'{u){v — u))wi = W1W2 = a.e. in Bi 

and, in particular, 

g{v) — g{u) — g'{u){v — u) = W2 = a.e. in {x G Bi : wi{x) > 0} . 

Since g is strictly convex, this implies that v — u = W2 = a.e. in {x G Bi : wi{x) > 0}. We 
have thus shown that v — u = a.e. in {x G Bi : wi{x) > 0} and, since v — u = W2 < a.e. 
in {x G Bi : wi{x) = 0}, we conclude that v — u < a.e. in Bi. 

(ii) As in ifTOl Lemma 1 1.3.4], we define 

/(^) := {-dv + (1 - ^)m, 0) + a /" g{^v + (1 - ^)u)(p dx , ^ G [0, 1] , 
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where is the unique positive eigenfunction of the linearized operator B/S? — TA + \g'{u) in 
Hf) r(®i) associated to the eigenvalue /Ui(n) = satisfying ||(/)||i = 1, see Lemma [Z9l Since g 
is convex, > in Bi, and -dv + {1 — d)u satisfies 

BA^{:dv + {l-^)u)-T/\{'dv + {l-^)u) + \((dg{v) + {l-^)g{u))<Q in Bi , 

we conclude 

/(^)<0, ^€[0,1]. (2.15) 

As 

/'(^) = (^; _ y, 0) + A f g'i'dv + (1 - ^)u){v - u)(t) dx 

and 

/"(^) = A / g'X'dv + (1 - ^)u){v - ufcj) dx , 



the assumption yUi(M) = guarantees that /'(O) = while the convexity oi g and the positivity 
of imply that /"(O) > 0. In addition, recalling that /(O) = 0, we deduce from (12.151 ) that 
/"(O) < 0. Therefore, /"(O) = and the strict convexity of g and the positivity of (j) in Bi 
entail v = u. D 

In order to study more precisely the behavior of solutions to S^ as the parameter A varies, we 
now derive several estimates. 

Lemma 2.11. There is Ci > such that 

WAh^ + ll^*llc3/2(Bi) + -^ / 9{u{x)) dx < Ci (2.16) 

whenever A G [0, A*] and u G S^. 

Proof. According to Lemma |Z2] there are mi > and 0i G C^{%i) satisfying Acjyi = mi0i 
and 

01 >0 in Bi, ||0i||i = l. (2.17) 

Multiplying (12.11 ) by (pi and integrating over Bi give 



—mi / (f)iudx = X g{u)(j)idx . 

Since u > — 1 in Bi, we deduce from (12.171) that 

< A / giu)(f)idx <mi . (2.18) 

Next, recall that Lemma 1277] ensures that the function u : [0, 1] — t- M defined by u(|x|) = u{x) 
for X G Bi is non-decreasing. This readily implies that g(u) is non-increasing and, thanks 
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todnn), 



0<A/ g{u{x)) dx = X\dMi\ giuir)^-' dr 



1 



/•3/4 /•! 

<A|9Bi|/ g {u{r))r'^~^ dr + X\dMi\ ^(u(r - l/4))/~Mr 
Jo J3/4 

/•3/4 1-3/A 

<A|9Bi|/ ^(u(r))r^-^dr + 2^-^A|9Bi| / ^(u(r))/-^ dr 

io Jl/2 

1-3/4 

< Xil + 2'^~^)\dMi\ g{uir)y-^dr 

Jo 

2'^A f 
< —. — / (t>i{x)g{u{x)) dx 

minB3/^ 01 iB3/4 

2'^mi 



mingg/^ 01 



We have thus proved that 

A||(7H||i<Ci. (2.19) 

It next follows from (I2.1I) . (I2.19I ). and the non-negativity of g and 1 + n that 

c||M||_ff2 < {u,u) = —A / g{u)udx < A||5'(n)||i < Ci . (2.20) 

J Ml 

Finally, if d = 1, the embedding of H'^(Ei) in C^/^(]Bi) completes the proof in this case. If 
c? = 2, we note that u solves BA'^u = TAu — \g{u) in Bi subject to homogeneous Dirichlet 
boundary conditions with ||TAn — Xg(u)\\i < C hy (12.191 ) and (|2.20l) . Hence, in this case the 
assertion follows from a version of the Brezis-Merle inequality [|4l (see Lemma IaTI) and the 
embedding of W^g^(Bi) in C^/^(Bi) for q large enough. D 

Restricting our attention to u E S^ with a non-negative /ii(u), the previous estimates can be 
improved in the following way. 

Lemma 2.12. There is C2 > such that 

/dx 

whenever A G [0, A^,] and u G S^ with /ii(n) > 0. 

Proof. We infer from (12.11 ) and the assumption yUi(n) > that 

—A / g{u)udx= {u,u) > —A / g'{u)u'^ dx 



A FOURTH-ORDER MEMS MODEL WITH CLAMPED BOUNDARY CONDITIONS 17 

and thus 

^?^ dx < . (2.22) 

Observing that ?,z'^ + z > 1/4 for 2 G (-1, -1/2), we deduce from (12.221) that 



/ 



1 dx< / '^:^/^'-;l"^ dx + / '^:^-^/^!l"^ dx 



^a+«)3 -4^ (l+«)3 ' y^^ (l+«^3 



, If-i _i/2-)(ti)(3n^ + n) 

< 8 li +4 ^ '• 7^^^ ^\, '- dx < 

(l + n)3 



Finally, (12.11 ) and Holder's inequality give 



\ 2/3 



0<{u,u) = -X ^(n)ndx<A,|Bi|^/M / -^ dx 

and (12.211) follows from the previous two inequalities and the finiteness of A^,. D 

2.3. Maximal stationary solutions. We first recall the existence of maximal solutions to (12.11 )- 
(IZ2I) established in [[M Theorem 1.5]. 

Proposition 2.13. (i) For any A G (0, A*), the set S^ is non-empty and contains a unique 
maximal element u\ in the sense that u < u\ for all radially symmetric classical subsolutions 
u to (|2.1I) - (I2.2I) with parameter A. In addition, for each x G Bi, the function A H- u\{x) is 
decreasing in (0, A*). 

(ii) There is no radially symmetric classical solution to (I2.1l) - (l2.2l) /c>r A > A^,. 

We supplement Proposition l2. 131 with continuity properties of A 1 — v u\. 

Lemma 2.14. The map A 1 — )■ u\ is continuous from [0, A*) to (7^(Bi) with uq = 0. In addition, 
the map A 1 — )■ jiiiux) belongs to C([0, A*)). 

Proof Fix A G [0, A^,) and let (Afc)fc>i be a sequence in [0, A*) such that A^ — )■ A as A; — )• 00. 
Then there is r/ G (0,1) such that 

A < -qX^ and A^ < -qX^, < X^, , k > 1 . 

Proposition 12. 131 ensures that ux^^ > Ur,x^ for all A; > 1, so that {uxf?)k>i ranges in a compact 
subset of (—1, 0]. Therefore, {g{uxj)k>i is bounded in Loo(Bi) and classical regularity results 
entail that {ux^)k>i is bounded in W^{Mi) for all q G (l,oo), see DUl Theorem 2.20] for 
instance. The compactness of Sobolev's embedding then implies that a subsequence of {uxf,)k>i 
(not relabeled) converges weakly in H'^(Mi) and strongly in C^(Mi) to a function u G Hf-, ^(Bi), 
which is a strong solution to (|2.1| )- (|2.2| ) and satisfies u > u^x, > —1 in Bi. Since g is smooth 
in (—1, 00), there is a > such that g{u) belongs to C^'^°'(Bi) and a further use of classical 
elliptic regularity results guarantees that u actually belongs to S,^, see flV, Theorem 2.19] for 
instance. 
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Consider now a radially symmetric classical subsolution a G C^(]Bi) fl C^(]Bi) to (I2.1I) - (I2.2I) 
with a > — 1 in Bi. For i) G (0, 1) and k > 1, we infer from the convexity of g and the 
properties Afc < 77 A* and g{urjx,) > 1 in Bi that 

BA^ ((1 - ^)a + ^u^aJ - TA ((1 - i})a + ^u^aJ + A,^ ((1 - ^)a + i}u^x,) 
<{l-t)) (5AV - TAa + Xkg{(r)) + ^ {BA\x, - TAu^x. + \kg{u,x.)) 
<(1 - ^){Xk - X)gia) + ^(Afc - vK)giu^x,) 
<\\k - A|||5f(cr)||oo -'i9(r/A* - A^) . 
Since A^ -> A as /c -> 00 and A < ?7A*, there is k^ > 1 large enough such that, for all k > k^, 

\\k - MM^)\\oo - HvK - \k) < , 

and hence (1 — ??)cr + ?9n^A, is a subsolution to (|2.1|) - (|2.2I) with parameter A^. Therefore, owing 
to the maximality property of ma^ , 

(1 - i9)a + i^Ur,x, < ux,, in Bi 

for all k > k^. We first let A; — ;■ 00 and then ^ — ?> in the above inequality to conclude that 
cr < n in Bi. In other words, n is a maximal solution to (|2.1I) - (|2.2I) and thus u = ux- 

Owing to the definition (I2.11|) . the continuity of A 1 — )■ Hi{ux) in [0, A*) readily follows from 
that of A I — > Ux which we have just established. D 

The next proposition entails that the maximal solutions are exactly the linearly stable solu- 
tions. 

Proposition 2.15. Let A G [0, A*). Then fii{ux) > 0, and ifu G S^ satisfies fii{u) > 0, then 

u = Ux. 

Proof. Due to the monotonicity and negativity of g' and the monotonicity of A 1 — > ux stated in 
Proposition |2?T3l it readily follows from (12.111 ) that 

fJ'iiux,) > /ii(nAj for < Ai < A2 < A=, . 

Introducing 

Xst := sup {A G [0, A*) : fii{ux) > 0} , 
we assume for contradiction that Xgf < K- Lemma [2. 141 then ensures that fii{ux^t) = 0. Now, 
given A G (A^i, A*), we deduce from (12.11) that 

BA\x - TAux + Xstgiux) = {Xst - X)g{ux) < in Bi . 

Applying Lemma r2 . 1 Ol (ii) . we conclude that u \ = ux^t and a contradiction. Therefore, Xst = K. 
Finally, considering u E S^ such that ij,i{u) > 0, Lemma [2.1 01 (i) implies ux < u while the 
maximal property of ux guarantees u < ux. Therefore, u = ux. □ 

We now show that the maximal arc-connected set ^0 defined in (12.61) coincides with the 
branch of maximal solutions (A, ua), A G (0, A*). 
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Corollary 2.16. There holds Aq = A* and V{X) = ux for each A G (0, A*). Moreover, 
\imx-^x, fii{ux) = 0. 

Proof. Since V{\) is for each A G (0, Aq) a classical radially symmetric solution to (I2.1I )- (I2.2I ). 
we clearly have Aq < A*. We now claim that 

/ii(V(A)) > for A G [0, Ao) . (2.23) 

Indeed, the continuity of V entails that A i-)- ^i{V{X)) is continuous on [0, Aq) with 

/ii(V(0)) = /ii(0) = mi>0 

with mi introduced in Lemma IZ2l Clearly, Lemma |2!9] (ii) implies fii{V{\)) > for each 
A G [0, Ao). Consequently, ux = V{X) for each A G [0, Aq) by Proposition l2.15[ 

Now, suppose for contradiction that Aq < A*. Then V^(A) — )■ uaq in C^{Mi) as A — )■ Aq 
according to Lemma 12.141 Owing to the maximality of Aq, this implies that F„(Ao,mao) is 
not boundedly invertible which contradicts Proposition 12. 151 since /ii(nAo) > 0. In particular, 
limA^A. l^i{ux) = due to the maximality A^,. D 

We finally investigate the behavior of ma as A — ;■ A^,. 
Proposition 2.17. For x G li, define 

Ux ix) := inf Ux(x) = lim Ux(x) G [—1,01 . (2.24) 

Ae(0,A.) A->A. 

Then ux, G S^* and Hi{ux,) = 0. Moreover, any u G S^* satisfies u < ux^ in Bi and if, in 
addition, fii{u) > 0, then u = Ux,. 

Proof. The proof is similar to that of [|TOl Theorem 11.4.1]. Indeed, the fact that ma, is well- 
defined is a simple consequence of Proposition |27T3](i). Thanks to Proposition 12. 15[ we are in 
a position to apply Lemma 12.121 and conclude that (nA)A is bounded in iJ^(Bi) while {g{ux))x 
is bounded in L3/2(Bi). Consequently, the sequence (nA)A is bounded in 14^3^2 (Bi) by classical 
elliptic regularity, see [ITTl Theorem 2.20] for instance, so that ux, G C^(Bi) due to the contin- 
uous embedding of Wy2(^i) i^i C^(IBi). If the minimum of ux, in Bi would be equal to —1, 
then Ma, (0) = — 1 and Vma. (0) = according to Lemma 12771 These properties entail that there 
is C > such that 

1 C - 

1 + Ux,[X) \x\ 

which contradicts the boundedness of {g{ux))x in L3/2(Bi). Therefore, ua, > —1 in Bi which, 
together with the above mentioned estimates and classical elliptic regularity, implies ux, G S^* . 
That ^i(ux,) = follows from Corollary 12.161 
Finally, take u G S^* . Then, for all A G (0, A*), 

BA'^u - TAu + Xg{u) = (A - X^)g{u) < in Bi , 

and thus n < ua in Bi by Proposition 12.131 Letting A — )■ A* gives u < ux^ ■ The uniqueness 
assertion is a consequence of Lemma 12.101 and the maximality of ma, • Q 
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2.4. Continuation of maximal solutions. We shall now completely characterize the radially 
symmetric solutions to (I2.1I) - (I2.2I) near (A*, uaJ. That is, we show that the curve A defined 
in (12.71 ) non-trivially extends the curve {(A, ux) : A G (0, A*)} of maximal solutions established 
in Subsection |2.3[ the latter coinciding with the curve ^o from (12.61 ) as proven in Corollary 12.161 
Moreover, all radially symmetric solutions to (|2.11 i- (|2.21 i near (A^,, ux,) lie on this curve, which 
in particular excludes any bifurcation phenomenon near this point. 

Indeed, let us first note that (A*, uxj cannot be the limit of (A(s), U{s)) as s — )■ oo owing to 
Proposition l2.17l and Theorem 12. 5 1 (iii). Thus, there is a minimal s* G [1, oo) such that 

(A(s,),t/(s,)) = (A„maJ and Ao = (A, t/)((0, s,)) . (2.25) 

Next recall from Proposition [OT] that ux, G S,^* with fii{ux,) = and it thus readily follows 
from Lemma |Z9] that the kernel of the partial Frechet derivative -Fm(A*, ux,) of the function F 
defined in (12.51 ) is spanned by a positive function 0* G Cfj ^(Bi). 

Now, a precise description of the behavior of A near (A*,maJ can be obtained from [[81 
Theorem 3.2] and is stated in the following theorem. 

Theorem 2.18. There are 5 > 0, e > 0, and an injective and continuous function Qfrom (—5, 5) 
onto (s* — 5, s* + e) with the following properties: 

(i) C(0) = ■%; 

(ii) (A, U) o C is a real analytic function on (—(5, 5); 
(iii) all solutions (A,m) to (|2.1| )- (|2.2| ) near the point (A*,maJ = (A(s^,), f/(s*)) lie on the 

curve {(A, U) o ((cr) : |cr| < 5}; 
(iv) (A o C)'(O) = and {U o C)'(O) = 0*. 

Proof. Since ux, G S^* we have ux, G O with O defined in (12.41 ) and thus the function F is an- 
alytic near (A*, maJ- As the kernel of F„(A^,, maJ is one-dimensional, codim(rg(Fu(A*, maJ)) 
equals 1 since -F„(A^,, maJ is a Fredholm operator of index 0, see e.g. [|22] Theorem 4.25]. We 
now claim that 

Fx{K,ux,) = A-^g{ux,) ^ rg(F„(A*, uaJ) . 
Indeed, if otherwise there is v? G C\) ^(IBi) with 

A~^giux,) = ip + KA~^g'{ux,)ip ■ 
But then ip G i/^^(li) satisfies 

BAV - TA(^ + A,^'(nAj<^ = g{ux,) in Bi , 
and testing this equation with 0* > yields the contradiction 

0= / [i?A^0^ - TA0* + A*5f'(MAj0*] "y^dx = / 5'(maJ 0* dx > . 

Therefore, we are in a position to apply [8, Theorem 3.2] and obtain in combination with The- 
orem |23]the assertion. D 

Actually, the curve A bends down at (A*, ux,)'- 
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Corollary 2.19. There holds (A o C)"(0) < 0. 
Proof. Twice differentiation of the equality 

= F(A o C(a), U o C(a)) = U o C(a) + A o ao^)A-^g{U o C(a)) 
with respect to a at cr = and Theorem 12.181 yield 

(A o C)"(0)A-i^KJ + KA~'g"{u^M = -F^iK, uxJiU o C)"(0) , 
that is, 

(A o C)"(0)^(«aJ + A./(«aJ0^ = -A{U o C)"(0) - X.g'{u,J{U o C)"(0) . 
Testing this last equation with 0* > and using the convexity of g imply 

(AoC)"(0) [ giux,)(l)*dx = -X, [ g"iuxMdx<0, 



whence (A o ()"(0) < 0. D 

2.5. End point. The following theorem now completes the picture of the curve A defined 
in (|2.7I) . It characterizes the limit of (A(s), U{s)) as s — )■ oo and shows that for each A G (0, A*) 
there are at least two steady-states. 

Theorem 2.20. (i) As s ^ oo, 

{A{s),U{s))^{0,uj) in Rx (C(ii) n C^(ii \ Ip)) 

for each p G (0, 1), where uj G C^(IBi \ {0}) fl C^(]Bi) solves the equation 

BA'^uj{x) - TAlu{x) = for x G Bi \ {0} , 

with boundary conditions 

uj{x) = dyUj{x) = for X G dE>i , 
a;(0) + 1 = Vcj(O) = , 

and satisfies u{x) > —IforO < \x\ < 1. 

(ii) For each A G (0, A*) there are at least two values < si < s* < S2 with A{sj) = A, 
U{sj) G S^forj = 1, 2, and [/(sa) < U{si) in Bi with [/(sa) 7^ U{si). 

Note that Theorem 12.201 allows one to compute the end point tu explicitly in terms of the 
modified Bessel functions of the first and second kinds for d = 2 and, respectively, in terms of 
the exponential function for d = 1 (cf. Figured])- 

To prove Theorem 12.201 we first need the following result relating the minimum of a function 
w to the integrability of g(w). 
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Lemma 2.21. Let w be a radially symmetric function in H^^.{Mi) such that g{w) G Li(Bi), 
and let the profile xvofw defined by tt)(|a;|) := w{x) for x E Mi be a non-decreasing function 
on [0, 1]. Then there is C3 > such that 

/d(l + wm < C3 (1 + Ml,) \\g{w)\U , (2.26) 

where d = 1,2 and 

Proof. Since ro is non-decreasing on [0, 1], we deduce from the Cauchy-Schwarz inequality 
that, for r e [0,1], 



< r'^-^a,tt)(r) = ds (s^"^9,tt)(s)) ds 
Jo 



d 



y^d I n- 







1 



jds (/-l9,tt)(s)) 



2 



1/2 



s'^-^ ds 



^d/2 

< ||A7/>|L . 

Then < 9rtn(r) < C|| Au7||2r*^^~'^^/^ and integrating once more with respect to r gives 

ra(r) < tt)(0) + C\\^w\\2 r^^"^^'^ , r G [0, 1] . 
Consequently, 

g{w{x))dx> I ^(m;(0) + C||Au;||2 Ixl^^-'^)/^^ dx . 

Setting ^ := (1 + w(0))^/'^'*~'^) and restricting the integral on the right-hand side of the above 
inequality to Bi \ Bg, we obtain 

q{w{x))dx > / 7T dx 

~ M\B, (1 + w{0) + C|| Au;||2 |a;|(4-rf)/2)' 

, 1 dx 



AB, (1 + C||A«;||2)' Ix]^^-"^ 

whence (IZ261) . D 

Lemma 2.22. Le? (A„)„>i Z?e a sequence of real numbers in [0, A*] and {vn)n>i be such that 
Vn G S^" for each n>l. If 

lim niinw„ = —1 , (2.27) 



r 

lim min '^i — — 1 
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then there are a subsequence of (A„, f „,)«>! (not relabeled) and uj G C^(IBi \ {0}) fl Cr(IBi) 
such that Lo solves 

BA^uj{x) - TAuj{x) = for x G Bi \ {0} , (2.28) 

uj{x) = d^uj{x) = for X G 9Bi , (2.29) 



and 



lim A„ = , (2.30) 



n—^oo 



lim \\vn - w||ci+cfB ^ == , a G [0, 1/2) , (2.31) 

w„ ^ w in H^{Mi) , (2.32) 

a;(0) = -lanJ uj{x) > -1 for x G li \ {0} . (2.33) 

Proof. Step 1: Compactness. By Lemma [2.11l and the finiteness of A*, {vn)n is bounded in 

if^(Bi) n C'^/^(Bi) and (A„)„ is bounded in [0, A*]. The compactness of the embedding of 
C'^/^(Bi) in C^"''"(Bi) for a G [0, 1/2) guarantees that, after possibly extracting a subsequence, 
we may assume that there are A^o G [0, A*] and a; G iJ^(Bi) n C^+°(Bi) for a G [0, 1/2) such 
that A„ — )■ Aoo as n — 7- oo and (12.311) and (12.321) hold true. Combining (12.271 ) and (12.311 ) readily 
gives 

uj{0) = -1 , Va;(0) = , and u = d^oo = on 9Bi . (2.34) 

In addition, by Lemma [2/71 f„ is radially symmetric with a non-decreasing profile D„ defined 
by Dn(|a^|) := Vn{x) for x G Bi. Consequently, the function cu enjoys the same properties by 
(12.311) and its profile D, defined by t>{\x\) := uj{x) for x G Bi, is a non-decreasing function on 
[0, 1]. Therefore, it follows from this property and (|2.34) that there is a G [0, 1) such that 

u}{x) = -1 for xeMa and u{x) > -1 for x G Bi \ !„ . (2.35) 

In addition, if a > 0, then 

d^uj{x) = for X G dMa . (2.36) 

Step 2: Identification of Aoo- To this end, we apply Lemma [2311 and use Lemma [TTTI to 
obtain 

XMl + v^m < Cs (1 + l|t;„||^0 A„||^K)||i < CiCsil + C^) , (2.37) 

which also reads 

^" (^TTTvfml^) - ^ '^ "^"^ ''''^ A„|ln(l + i;„(0))|<C if d = 2. 

Letting n — )■ oo in the above inequality readily gives Aoo = by (|2.27l) . whence (12.301) . 

Next, fix g E (a, 1). We infer from (12.311 ) and (12.351 ) that {g{vn))n>i is bounded in Loo(IBi \ 
Mg), SO that {Xng{vn})n>i convcrgcs to zero in Loo(Bi \ B^) as n — )■ oo by (12.301) . Classical 
elliptic regularity estimates then allow us to pass to the limit as ra — )■ oo and conclude that 
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u G C^{Mi \ Bg) satisfies B^^uj{x) - TAuj{x) = for all x G Bi \ B^ and (IZ291) . Since g is 
arbitrary in (a, 1), we have shown 

BA^uj{x) - TAuj{x) = for x G Bi \ B„ . (2.38) 

Step 3: Identification of a. The final step is to prove that the yet unknown number a is equal 
to zero. 

Let n > 1. According to (12.101 ). there is r„ G (0,1) such that 

At;„ < in Bi \ B^„ and Av^ > in B^„ . (2.39) 

The boundary conditions for Vn and Lemma [2.11| then imply 

|At>„(x)| dx = — Avn{x) dx = — Af„(x) dx + Avn{x) dx 

< v^Bj ||At;„||2 < Civ^Bj . (2.40) 

Assume for contradiction that there is a subsequence (r„^)fc of (r„)„ such that r^^. — )■ as 
A; — 7- oo. It readily follows from (12.401 ) that {Avnf,)k converges to zero in Li(Bi \ B^) for all 
Q G (0, 1). Recalling (I2.32L we deduce that Acu = almost everywhere in Bi which, together 
with (12.291 ). implies a; = in Bi and contradicts (12.351 ). Therefore, there is r^ > such that 

rn>n>0, n > 1 . (2.41) 

Now, forn > 1 and x G Bi, we set 0„(|x|) = Vn{x), rD„(|x|) := Avn{x), and define 

an := sup Ir G (0, r^) : 0„(r) < -1 + a/A^| 

if the set is non-empty, and (t„ = otherwise. Since (t„ G [0, r^], we may assume, after possibly 
extracting a subsequence, that 

Urn an = a e [0,n] . (2.42) 

Assume a > for contradiction. The definition of an, (12.311) . (12.391) . and (12.411) then ensure 
that 

uj{x) = -1 for X G B^ and Atv(x) > for x G B<,„ . (2.43) 

Since 1 + Vn{x) < a/A^ for x G Bo-„ and Vn G 5^", we find 

BA^n - TAvn = -\ng{vn) < "1 in B,„ . 
Consequently, 

dr (r'^-i {BdrWnir) - TdrMr))) < -r"-^ , r G (0, an) , 

and, because r'^~'^9rtt)„(r) and r'^~^9rti„(r) both vanish as r — )■ 0, a first integration gives 

dr {BWn{r) - Tt)„(r)) < -^ , r G (0, (T„) . 
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We next integrate the above differential inequality over (r, (T„) to obtain 



(jH r"^ 



5m„(r) > BxOnM + T (0„(r) - 0„(o-„)) + ^ - ^ 



^2^d—l j,d+l 



Due to rD„(a„) > by (12.431) . we find 

Bdr (/-^9,0„(r)) = 5r"-^tt)„(r) > Tr"-^ (D„(r) - 0„(a„)) , ^^ 2rf , 

whence, after integrating once more, 

Br''-'drX>^{r) >T j^^ s''~' (o„(s) - 0„(c7„)) d^ + ^ - ^^|^ ' '^ ^ (0, a„) . 

Now, fix r G (0, a). Owing to (12.311 ) and (12.421) , we may pass to the limit as n — )■ oo in the 
above inequality and deduce 

Br^-%Mr) > T [ .-' („(.) -„(.)) d. H- ^ - ^£L^^ . 

As r e (0, cr), it follows from (12.431) that 9r.t)(r) = and t)(s) = -1 = D(cr) for s e (0, r), so 
that we end up with 

= Br'-'dMr) > ^^^^^ [{d + 2)a' - dr'] > , 

and thus a contradiction. We have thus shown that 

lim a„ = . (2.44) 

n—^oo 

Let then r E (0, 1) be arbitrary. Owing to (I2.44I ). there is Nr > I large enough such that 
an G (0, r) for n > Nr. Recalling the definition of an, this means that, for n > N^ and 
x e Bi \ Mr, we have 1 + w„(x) > a/A^ and thus \ng{vn) < 1 in Bi \ B^. Since (f„)n>Af,, is 
bounded in if|,(Bi) by Lemma [UTI and Vn G Sr", we conclude that {vn)n>Nr is bounded in 
if|,(Bi \ Br). The convergence (12.311 ) then entails 

oj e HfyiMi \ Mr) for all r G (0, 1) . 

Assume now a > for contradiction. Since a; = — 1 in B^, the just established regularity of u 
leads us to 

= 1 + 0(a) = drti{a) = 9>(a) = d^t){a) . (2.45) 

Multiplying (12.381) by u, integrating over Bi \ !„ and using (12.291 ) and (12.451) give 






This implies that uj is constant in Bi \Ba and contradicts (12.291 ) and (12.451 ). Consequently, a = 
and the proof is complete. D 

To finish off the proof of Theorem 12. 201 it just remains to summarize our previous findings. 
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Proof of Theorem^jUi Recall from Theorem |231 (1X25]) . Corollary llHl and <^M that A = 
{(A(s), U{s)) : s > 0} is a continuous curve with 

lim(A(s), U{s)) = (0, 0) , (A(s,), U{s,)) = (A„ n^J , 

s—>-0 

and ^0 = {(A, Ma) : A G (0, A^,)}. Part (i) of Theorem 12 .20 1 now follows from Lemma 12.221 
and Theorem 12.51 (iii). 

Since A(s) — > for s -H- oo as just shown, we find for each A G (0, A*) numbers < si < 
s* < S2 depending on A such that A(si) = A(s2) = A. Since si < s*, we have U{si) G ^o and 
thus f/(si) > U{s2) by Proposition [Ol Corollary [IHl and (IZ25T) . Moreover, U{si) ^ U{s2) 
since no bifurcation can occur along the curve ^o due to the implicit function theorem. This 
proves part (ii) of Theorem 12.201 D 

Theorem 11.11 and Remark 11.21 are now consequences of Theorem 12. 5[ Theorem 12.181 and 
Theorem 12.201 

3. Well-posedness in general domains 

We shall now focus on the well-posedness of the dynamic problem. Let us recall that Vt is an 
arbitrary smooth domain in W^ with d = 1,2. 

3.1. Well-posedness for the hyperbolic problem. In this subsection, we prove Theorem [Ol 
To lighten the notation, we agree upon setting 7 = 1 in the following. We first reformulate 
(ll.ll) - (ll.3|) as a first-order Cauchy problem and use well-known results on cosine functions 
(see e.g. [HI Section 5.5 & Section 5.6] for details). Let us note that the self-adjoint operator 
—A = —BA^ + TA with domain Hj^iVt) generates an analytic semigroup on ^2(^2) with 
spectrum contained in [Re 2; < 0], its inverse A'^ is a compact linear operator on ^2(^2), and 
the square root of A is well-defined. Noticing that A is associated with the continuous coercive 
form 



{u,v) 



I {BAtiAv + TVn ■ Vv) dx , u,v e Hl{n) 
Jq 



the domain of the square root of A is (up to equivalent norms) equal to Hj^iVt). Consequently, 
the matrix operator 

'0 -T 



^-•A 1 

with domain D{K) := Hj^iVt) x Hj^iVt) generates a strongly continuous group on the Hilbert 
space HI := H'j-,{yt) x L2(f2). Writing uq = (n°, u^), u = (n, dtu), and 

/(u)=f_l,^l with g{u):=l/{l + u)\ 



we may reformulate (|1.1I )- (|1.3I ) as a Cauchy problem 

u + Au = A/(u), t>0, u(0) = uo (3.1) 
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in H with li indicating the time derivative. Now, defining for k E (0, 1) the open subset S{k) of 

Hby 

S{k) := {u G Hl{n) : u> -l + nm^jx L2(fi) , 

the function f : S{k,) —^ Mis uniformly Lipschitz continuous. A classical argument then entails 
the following proposition. 

Proposition 3.1. For each uq G S{k), the Cauchy problem (13.11) has a unique maximal mild 
solution vi = {u.dtu) G C{[0,Tm),M.) for some maximal time of existence Tm = Tm(uo) G 
(0, oo]. Ifvm < C)0, then 

liminf ( mm u(t)) = — 1 , (3.2) 

or 

lim sup \\{u{t),dtu{t))\\M = oo . (3.3) 

To obtain more regularity on the mild solution u, let us consider an initial value in the domain 
of the generator —A, that is, let uq G (if^(f2) x if|,(f2)) fl S{k). Then, since / is Lipschitz 
continuous, it follows as in the proof of [|20l Theorem 6.1.6] that u : [0, Tm) — )■ H is Lipschitz 
continuous and whence differentiable almost everywhere with respect to time. Consequently, 
we obtain (see also [|20l Corollary 4.2.1 1]): 

Corollary 3.2. If uq G (^Hf){Q) x _ff|,(r2)) n S{k), then the mild solution u is actually a strong 
solution to (13.11 ). That is, u is differentiable almost everywhere in time with u G -^i(0, r; EI)/or 
each T G (0, r^) and 

n{t) = -Au{t) + /(u(t)) 
in Mfor almost every t G [0, Tm). 

As a consequence, since u = {u,dtu),we deduce under the assumption of Corollary I3.2l that. 
for each r G (0, Tm), 

d^u G C{[0,Tm),Hl~"'{Q)) , d^+'u G L,{0,t;HI~'\Q)) , 

for /c = 0, 1 and 

(BA^ - TA)u = -d'fu - dtu - A(l + m)-^ (3.4) 

Since the right-hand side of (13.41) belongs to lvi(0, r; L2(fi)), we deduce u G i^i(0, r; HjjiVL)). 
Now, testing (Il4l) by dtu G C([0, r^), L2(ri)) results in 

-^( f IdM'^ dx+ [ (51 AmP + Tl VmP) dx - 2A /" -^ dx] = - [ l^iMp dx (3.5) 
2 dt Vif7 Jn Jn^ + u J Jn 

almost everywhere in [0, r^). Assume now that r^ < oo and that (13.21) does not occur. Then 
(1 + u)~^ G Loo((0,rm) X r2) so that (|3.3I) cannot occur as well, whence a contradiction. 
Consequently, r,„ < oo implies (13.21) . 
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To finish off tlie proof of Tlieorem 11.31 it remains to sliow tliat the solution exists globally 
in time for small A and small initial values. Recall that H'^(p.) embeds continuously in Loo(f^) 
since d = 1,2 and let C4 > be such that 

Ml < c4B\\Av\\1 + T\\Vv\\l) , V e HUQ) . 

Then we can prove the following result on global existence: 

Corollary 3.3. For each k G (0, 1/2), there exists \i{k) > such that r^ = 00 provided that 

X < Ai(k) and uo G (i/|,(fi) x Hl{n)) n S{2k) with 

BWAux + nvuji + hX < ^-^—^ . 

C4 
Proof. Since n° > — 1 + 2k, we have 

To := sup{r G (0, r^) : u{t) > -1 + k , t G [0, r)} > 

and (1 + u{t))-^ < K"i for t G [0, Tq). From (l33T) . 

o\ ini 
B\\Au{t)\\l + T\\Vu{t)\\l < BWAu'^Wl + T||Vn°||^ + Hn^H^ + -^-^ (3.6) 

K 

for t G [0, To) and therefore 

Wumlo < (1 - 2Kf + ^^^ < (1 - «:)^ te [0, To) , 

K 

if A < Ai(k) with Ai(k) > sufficiently small. Consequently, To = r^ from which t^ = 00 
by Proposition [3TI D 

Note that u G Too(0, cx); H^{^)) and M(t) > -1 + /t for t > due to (1X61 ) and To = cx). 

Remark 3.4. If Q = Bi, f/ze rotational invariance of (11.11) anJ ?/ze uniqueness of solutions 
guarantee that u(t) is radially symmetric for each t G [0, Tm) provided that {vP, u^) is radially 
symmetric. 

Remark 3.5. The proof of Theorem 17.31 is the same if the clamped boundary conditions (11.21) 
are replaced by the pinned boundary conditions (|1.6I) . and we obtain a strong solution in this 
case as well. This improves the existence result for weak solutions in W2\ . 

3.2. Well-posedness for the parabolic problem. To prove Theorem [L4l we first note that 

g:{ue HliVt) : u>-l + KmVL}^ H'^iVt) , m ^ (1 + m)-^ 

is uniformly Lipschitz continuous and recall that for instance H'^iVt) ^^ Hj^ (fi). We then 
also recall that the operator —A = —{BA^ — TA) with domain Hj^iVt) generates an analytic 
semigroup {e~*^ : t > 0} on L2{Vt) with 

||e-*^||£(L.(n)) < Me-°* , t>0, 
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for some a > 0. Formulating (ll.7l) - (ll.9l) by means of the variation-of-constant formula 

u{t) = e-*^u° - A / e-(*-^)^ g{u{s)) ds , t > , 



^0 

the proof of Theorem II ^I can be performed by a classical fixed point argument. In particular, 
the exponential decay of the semigroup entails global existence for small values of A as stated 
in Theorem II .41 (111) (e.g. see [|9l Theorem 1.2 (i)] for details). As in Remark |X4l if fi = Bi, 
we easily see that u{t) is radially symmetric for each t E [0, r^) provided that n° is radially 
symmetric. 

4. Touchdown in the ball 

In this last section we return to the case fi = Bi and take advantage of the fact that a positive 
eigenfunction 0i > to the operator A is available, see liT6l and Lemma [Z2l We employ the 
eigenfunction method as e.g. in ^ to show the occurrence of a singularity in finite time as 
stated in Proposition [T3] and Proposition 1 1.6 1 

4. 1 . General initial conditions. 

Proof of Proposition \1.5\ Recall that Lemma f22\ ensures the existence of mi > and (pi E 
Cj) r(Bi) with (pi > in Bi, ||0i||i = 1, and A(pi = mi(pi. Let u be the maximal solution on 
[0, Tm) to (fTTT])-(fr3l) if 7 > or (fT77Tl-(fr9]l if 7 = corresponding to the initial value (m°, u^) 
and define, for t E [0, r^), 



N{t) := (pi{x)u{t,x)dx>- (/)i(x)dx = -l. (4.1) 

Assume that 

Ami 
A > ^ . (4.2) 

We multiply (11.11) by (pi, integrate over Bi, and use the properties of (pi, the convexity of g, and 
Jensen's inequality to obtain 

7'^ + ^ < - /" {BA^(Pi - TA(Pi)u dx-XgU (Pm dx\ 

< - miN - \g{N) . (4.3) 

Setting x(^) := miz + \g{z) for z E (—1, 00), we note that 

X is decreasing in {—l,z\) and increasing in {zx,oo) , 

where zx := (2A/mi)^/^ — 1. The choice A > 4mi/27 guarantees x{^) ^ xi^x) > for 
z E (-1, 00). We then infer from (l4Tt and (|43T ) that, for t E [0, t^), 

dN 

-^it)<-xizx) and - 1< iV(t) < iV(0) - x(^A)t , 
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if 7 = 0, respectively 







^W.e-- 




- 


- Xizx) 


and 


r T TV T 1 






-1 <N{t) <N{0)+-f^ 


^(0)..(.) 


(- 


- e^"^' 


if7>0 


Consequently, 






Tm < 


l + iV(0)+72 


v 


>..( 


^a) 


Y 


1 

X{zx) 



X{z\)t 



< oo . 



This completes the proof of Proposition 1 1.51 



n 



4.2. Radially symmetric initial data. Roughly speaking, the proof of Proposition [L6] pro- 
ceeds along the same lines of that of Proposition [T3] but takes advantage of the properties of the 
linearization of (12.11 ) for A* described in Lemma |Z9] and Proposition 12.171 

Proof of Proposition \1.6\ Fix A > A* and let u be the maximal solution on [0, Tm) to (11.11) - 
(11.31) if 7 > or (ll.7l) - (|1.9l ) if 7 = corresponding to the initial value (m°, u^). Recall that 
/^i(waJ = by Proposition 12. 171 and that there exists a corresponding positive eigenfunction 
0* G Cj) r(IBi) to the operator A + X^,g'{ux,) according to Lemma |Z9l which we normalize so 

that \\(j)*\\i = 1. For t G [0, r^), define 



M{t) := / u{t,x)(f),{x)dx> 



(f)^{x) dx 



-1 . 



(4.4) 



As in IIT2I Theorem 4.1], we multiply (11.11 ) by 0^,, integrate over Bi, and use the equation 
satisfied by ux, to obtain 



,d^M dM 



7 



dt2 



dt 



+ 



[U 



, {BA\ - TAu + Xgiu)) dx 

* [BA'^ux, — TAux, + Kg{uxS) dx 

ux,) {BA^<p,-TA<i),)dx 

0* {Xg{u) — X^:g{ux,)) dx 

0* [-Xg{u) + X^g{ux,) + X^g'{ux,){u - uaJ] da; . 
It follows from the convexity of g and Jensen's inequality that 

X*g{ux,) + X^g'{ux,){u - maJ < X^g{u) and / g{u)(f)^ dx > g{M) . 
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Therefore, owing to the positivity of 0* and A — A*, 

.d^M dM 



<-\x-\,)g{M). 

Since (A — \^)g{M) > 0, a first consequence of (14.51) is that, for t G [0, Tm), 

dM 



(4.5) 



dt 



■(t) < and M{t) < M(0) 



if 7 = or 



,t/,^dM 



dt V dt 

if 7 > 0. In both cases, for t e [0, Tm), 



(t)] < and M{t) < Af (0) + 7'^(0) (l - e"*/^' 



M(t) < i^o,7 := M(0) + 7^ 



dM 



dt 



(0) 



Recalling that g is decreasing, we deduce from (14.51 ) and (14.61 ) that 



.d^M dM 



7 



whence, for t e (0, r^), 



if 7 = and 



^^2 + ^ + (A - A,)<7(i^o,7) < 
M(t) < M(0) - (A - X,)g{Kor/)t 
dM 



M{t) < M(0) - (A - A.)(7(J^o,7)^ + 7M -JT (0) + (A - A*)^(i^o,7) U - e^^^^ 



dt 
if 7 > 0. We have thus shown that, for t G [0, r^), 

M{t) < /fi,^ - (A - A,)(7(iro,^)t 
where 



dM 



dt 



(0) + (A - K)g{Ko,^) 



i^i,^:=M(0)+72 

Recalling that M (t) > -1 for all t G [0, r„) by (|44l) . we end up with 
-1 < i^i,7 - (A - K)g{Ko,y)t , t G [0, r„ 
Consequently, 

^ 1 + i^l,7 

(A - A,)5f(Ko,7) 
as claimed in Proposition 1 1.61 



(4.6) 



n 
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Appendix A. A Brezis-Merle estimate 

We shall prove here in two space dimensions that solutions to the biharmonic equation with 
homogeneous Dirichlet boundary conditions and right hand sides in Li belong to W^ for any 
q E {1, oo) (a fact which is used in Lemma [2.11|) . This result is strongly reminiscent of the cel- 
ebrated Brezis-Merle inequality dU stating that solutions to the Laplace equation with homoge- 
neous Dirichlet boundary conditions and right hand sides in Li belong to Lg for any g G (1, oo). 
The proof of Lemma IaTI below is actually very similar to that of 1*41 Theorem 1] and is given 
merely for the sake of completeness. Let us point out that, because of the clamped boundary 
conditions (11.51) . the result cannot be deduced directly from BH Theorem 1] (in contrast to the 
case of pinned boundary conditions (11.61 )). 

Lemma A.l. Let d = 2 and f G L2(Bi), / ^ 0. Let w G H^{Ei) be the unique solution to 

A'^w = f in li , w = duw = on dMi . (A.l) 

There are ^o > and C5 > independent of f and w such that 

exp — — dx < 65 , 



where D^w is the Hessian matrix ofw. Furthermore, given g G [1, 00), there is CQ{q) > 
independent of f and w such that 

\\w\\w^<Ceiq)\\fh. 

Proof. A classical density argument allows us to assume that / G C^{Mi). Introducing the 
Green function G associated with the operator A^ subject to homogeneous Dirichlet boundary 
conditions in Bi, it follows from [[TTl Theorem 4.7] that there is i^o > such that 



\d^,d:,^G{x,y)\ < i^oln 2 + mm <^ 1, 

V \x-y\ J [ \x-y[ 



for (x, y) eMi X Bi with d{-, 9Bi) denoting the distance to 9Bi. From this we deduce 

|<9,,9. ,G(x, y) I < Ko In (2 + -^—-] , (x, y) G Bi x Bi . (A.2) 



\x-y\ 
The solution w to (lA.ll) can be written as 



w{x) = / G{x,y)f{y)dy , x G Bi , 

which further gives 

d,^d,^w{x)= f d^d,p{x,y)f{y)dy, x G Bi , 1 < 2,j < 2 . (A.3) 
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Now, let 19 G (0, 2/Ko). We argue as in the proof of [4, Theorem 1] and use (|A.21 ). (IA.3t . the 
convexity of 2; i— )■ e"^^, and Jensen's inequality to obtain 

/■ f^\d^A.w{x)\\ f / /■ , , ,1 /(y) \ 

</ / exp(tf|a..a,jG(rc,!/)|)jfMd!,dl 

JBiJMi II J 111 

Jbi ll/lli Jbi V k-z/l/ 



For (x, y) E Ml X Ml, we have a; G 182(1/) and thus 



i9i^o 



exp ^l^f;fWI d.</^; 1^) d.d,<C(*). (A.4) 

since — ^A'o > — 2. 

Consider next g G [1, 00). We infer from (IA.4I I that 

'^J"'llg — 11^ 111 / II #l|9 



5..5.,^||!<ll/ll?/ ' T.nr' dx 



^Bi IIJ 111 

<||/||? sup{.^e-^/-o} / exp f 'Yimi"^' ) ^" 

^>0 JBi \ -'^oil/Ill J 

<C(l/iro) 11/111 sup {z%-^/^°}, 

z>0 

which, together with the Poincare inequality, completes the proof. D 

Remark A.2. According to [[TTl Proposition 4.27], the estimate (IA.2I) is valid for an arbitrary 
smooth domain Vt of M? (with a constant depending on Vt) so that the validity of Lemma IA.il 
extends to arbitrary smooth domains ofM.^. 
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